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An unsteady model is developed for the particle size distribution in fluidized-bed reactors 
including fragmentation, abrasion, elutriation and the chemical reactions of wood gasifi¬ 
cation. Based on the assumption of constant conditions (gas composition, temperature, 
velocity) of the surrounding atmosphere, an analytical solution is developed for the dis¬ 
tribution of sizes belonging to the classes of mother and fine particles. It is found that for 
the typical feed sizes (minimum above 3 x 10~ 2 mm) and the usual maximum size of fine 
particles (2.4 x 1CT 3 mm), the behavior of fine particles is quasi-steady with respect to 
mother particles. The numerical solution of the quasi-steady formulation of particle pop¬ 
ulation balances is also coupled with a two-phase (bubble and emulsion), three-zone (bed, 
splash zone and freeboard) model for a bubbling fluidized-bed reactor, giving predictions of 
the producer gas composition in agreement with measurements for air gasification of 
wood. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Biomass is the only widespread source of renewable energy 
and carbon. Increased use of this fuel can contribute in 
reducing carbon dioxide emissions and enhance energy se¬ 
curity. Thermo-chemical conversion has attracted significant 
interest, in particular gasification, as homogeneous, gas phase 
combustion of syngas can be more efficient and easier to 
control than solid fuel combustion. Fluidized-bed gasifiers 
present the advantages of good temperature control, low- 
medium tar yields, toleration of variation in the feed quality/ 
size, partial load operation, easy start-up and shut-down and 
high conversion efficiency although the particulate content in 


the gas and the carbon loss with ash are higher than for fixed- 
bed reactors [1], As pointed out in Ref. [2], the implementation 
of commercial plant technology traditionally foresees 
comprehensive experimental investigation, progressing from 
laboratory scale test units to a pilot-scale plant, before a 
commercial demonstration plant can be built. In order to 
optimize the process, extensive experimentation is required 
for the various scales about the effects of operating and design 
parameters. This optimization procedure can be speeded up 
and become less expensive with the aid of mathematical 
models, on condition that the equation formulation includes 
all the important chemical and physical processes by 
describing their mutual interactions and the dependence on 
the process parameters. 
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Several transport models have appeared [3-12] for biomass 
gasification in bubbling fluidized beds, which couple the well 
known two-phase hydrodynamics with global kinetics, mass 
and, in a few cases, heat transfer. However, none of them 
takes into account the distribution of bed particles (particle 
population). This affects and is, in its turn, affected by the bed 
hydrodynamics, the extent of chemical reactions, attrition, 
elutriation and the heat and mass transfer rates. 

A great effort has been made in the formulation and so¬ 
lution of particle population balances essentially in the field of 
coal conversion [13-19] with the exception of Khan et al. 

[ 20 . 21 ] who consider a fluidized-bed biomass combustor, but a 
comprehensive unsteady model is still lacking. The models 
currently available are for a large part for steady-state condi¬ 
tions [13,14,16-18,20,21] using several simplifying assump¬ 
tions. The unsteady models are based on either an analytical 
solution [19] or a numerical solution by means of the orthog¬ 
onal collocation method [15] but attrition of the particles is 
never taken into account for char, although the treatment 
proposed for other materials, such as limestone 22], could be 
extended to the problem under study. 

Two main mechanisms of size reduction for carbonaceous 
particles in fluidized beds are reported [23]: abrasion, where 
particles of much smaller size break away from the original 
particle, so that particles become slightly smaller than the 
original ones, and fragmentation, where the breaking-away 
process gives rise to a number of particles of smaller size 
than the original ones. Biomass particles and, in a special way, 
char particles, generated from the pyrolysis stage, undergo 
size reduction not only as a consequence of the reactions but 
also owing to fragmentation and abrasion. The attrited fines 
can be easily elutriated away unreacted, reducing the process 
efficiency. For instance, it is pointed out [17] that if a 
combustor is optimally designed or operated, each 1 % effi¬ 
ciency improvement would mean a tremendous cost saving 
for a power plant. 

This study consists of two main parts. In the first part, a 
comprehensive unsteady model of particle population bal¬ 
ance in fluidized-bed gasifiers is proposed. The particle dis¬ 
tribution is divided into two size ranges: the class of fine 
particles contiguous to the class of mother particles. This 
approach is supported by a number of previous studies 

[17.18.20.21] and is justified by the fact that there are no sig¬ 
nificant amounts of intermediate-sized particles between the 
mother particles and the fines produced from the former by 
abrasion. Based on the assumption of constant conditions of 
the surrounding gas, analytical solutions are developed of the 
model equations for the two size ranges which are used for a 
parametric analysis with the scope of identifying the pro¬ 
cesses/parameters which play a controlling role. Also, the 
conditions are determined which permit a quasi-steady 
approximation for the behavior of fine particles with respect 
to the class of mother particles. In the second part of the 
study, using the analytical solution as a benchmark, a nu¬ 
merical solution is implemented of the quasi-steady version 
of the model. This is incorporated in a state-of-the-art model 
for an isothermal bubbling fluidized bed biomass gasifier 
schematized according to the two-phase theory of fluidiza¬ 
tion. The coupled particle population balance and the fluid¬ 
ized bed models are used to simulate gasification of wood with 


air at a laboratory scale obtaining good agreement between 
measured and predicted gas composition. 


2. Mathematical model 

The unsteady, one-dimensional mathematical model for the 
bubbling fluidized bed is formulated in accordance with the 
two phase theory of fluidization [24], with the existence of an 
emulsion and a bubble phase, and introducing three spatial 
zones corresponding to bed, splash zone and freeboard, which 
are isothermal. The gasifying agent is injected from the bot¬ 
tom while the biomass particles are fed at a certain bed height. 
The gas phase behaves according to the ideal gas law and is at 
a constant pressure. The emulsion phase, consisting of an 
interstitial gas phase at the conditions of minimum fluidiza¬ 
tion and a solid phase, is assumed to be perfectly mixed. The 
fuel particles do not contribute in the definition of the void 
fraction and solid mass of the emulsion phase, given that they 
are only a small percentage compared with the inert particles. 
However, their size distribution (and mass) is determined by 
means of a population balance. The bubble phase consists of 
gas only and is described by a plug-flow approximation. 

A plug-flow approximation is used for the splash and 
freeboard zones. These, in addition to the gas phase, foresee 
the presence of a solid phase due to char particle elutriation 
from the emulsion phase: the particle size is constant and 
equal to the average value over the size range of the elutriated 
particles. Also, given the small quantity of solid particles, a 
unit void fraction and a velocity of the solid coincident with 
that of the gas are assumed. 

2.1. Chemical reactions and moisture evaporation 

Biomass pyrolysis gives rise to a huge number of chemical 
compounds which, for engineering applications, are often 
lumped into three groups: permanent gases, tar and char. 
They result from both primary decomposition of the solid fuel 
and secondary reactions of vapor-phase organic products into 
low-molecular weight gases and char, as they are transported 
through the particle and the reaction environment. The ki¬ 
netics of biomass pyrolysis has been extensively investigated 
[25] but, as the characteristic times are shorter by several 
order of magnitude with respect to those of the heterogeneous 
reactions of char conversion [26], the process is often assumed 
to occur instantaneously. This assumption is retained here, 
approximating primary degradation of wood by means of a 
global reaction: 

Wood-^cChar + vcoCO + v C o 2 C0 2 + v H U 2 + RCH 4 CH 4 

(al) 

+ hh 2 oH 2 0 + utiTj 

where the gaseous species, water vapor, primary tar (TJ and 
char (Char) are formed with assigned stoichiometric co¬ 
efficients. The stoichiometric coefficients (v C o = 0.035, 
i'co 2 — 0.035, i'ch 4 — 0.015, R 1120 — 0.20, V 112 = 0.0015, vq — 0.215, 
i'ti = 0.50) are those corresponding to the yields (%, wt) ob¬ 
tained from the pyrolysis of wood carried out in fluidized beds 
at temperatures around 800 K [26]. Instead, the thermal 
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cracking of primary tar, to give gaseous species, water and 
refractory tar (T 2 ), occurs with a finite rate: 

Ti >v co^ ,< ~ > + ''co 2 *-'® 2 + + i'J H4 CH 4 + vJj 2 CH 2 + i*t 2 T 2 (a2) 

t a 2 = A a2 exp ^ Pti (!) 

again with assigned values for the stoichiometric coefficients. 
The selected tar cracking kinetics is that reported in Ref. [27]. 
The stoichiometric coefficients are derived from Ref. [28] with 
some modifications. It is assumed that a single fraction of 
primary tar is formed. Moreover the yield of refractory tar 
(5.8%) is derived from [29], so that the exceeding fraction with 
respect to the figure given in Ref. [28] (17.5%) is re-distributed 
over CH 4 , H 2 and H 2 0. This treatment produces the following 
values for the stoichiometric coefficients: v* co = 0.46, 
v*co2 — 0.08, f*cH4 — 0.157,1**1120 = 0.23, r*h 2 — 0.015, i*t 2 — 0.058. 

Moisture evaporation is a low temperature process and 
again, given its short characteristic times compared with 
those of biomass pyrolysis and char conversion, is assumed to 
occur instantaneously at the feeding point. 

Char particles experience the heterogeneous reactions of 
combustion and gasification by means of carbon dioxide and 
steam: 

CH.,0^ + T 0 2 (2 - 2y-/? + ") CO + (2y+ 1) C0 2 + |H 2 0 

(bl) 

CH l( O„ + CO 2 ^2CO + 0H 2 O+ (|-/3 )h 2 (b2) 

CH„0 (i + (l- / 3)H 2 0 fe 4c0+ (l —/J + !)h 2 (b3) 

The elemental composition of char consists of 95% C, 2% H 
and 3% O (a = 0.2526, fi = 0.0237) [30,31] and the char/oxygen 
reaction produces a ratio CO/C0 2 = 0.5 (y = 0.8013). 

Some considerations should be made in relation to the 
selection of the conversion model for the char particles. The 
approximation of reacting particle (small Thiele numbers, 
regime I) is never achieved for practical conversion conditions 
[26], In this case, based on the evaluation of the Thiele num¬ 
ber, it is likely that combustion mainly takes place at the 
surface (regime III, shrinking particle) while gasification is 
confined to a surface layer (regime II, variations in both size 
and density of the particle) [26]. When dealing with reactor 
models, the description of single particle effects is always 
based on simplified treatments which, apart from a very few 
exceptions (for instance a combination of regime I for gasifi¬ 
cation and regime III for combustion [32]), uses the shrinking 
particle model (regime III) (for instance 
[4,7-9,20,21,30,31,33,34]) in some way assuming that the sur¬ 
face layer, where gasification reactions take place, is much 
smaller than the particle size so that all the reactions can be 
lumped at the surface. On the other hand, it should also be 
noticed that population balances, based on particle size dis¬ 
tribution, cannot handle the simultaneous variations in space 
and time of the particle density. Therefore, in this study, 
single particle effects are again described in accordance with 
the shrinking particle schematization (there is no surrounding 


ash layer, justified by the small ash contents of biomass fuels). 
To account for the simultaneous effects of diffusion through 
the gas film, adjacent the particle, and intrinsic chemical ki¬ 
netics, an effective reaction rate, r J( is introduced assuming a 
linear dependence on the oxidizing/gasifying species con¬ 
centration [30,31]: 



5j = y,l,(l-/3); i 


= Aj exp 



j = bl, b2, b3; 


0 2 , C0 2 , H 2 0 


(2 -3) 


where as already done in Refs. [7-10,35], the global mass 
transfer coefficient is described by the Ranz-Marshall corre¬ 
lation [36]: 




Sh ; = 2 + 0.6Sc° 33 Re 05 , 

p B u m fd c 

Re = —- 


(4-7) 


The kinetic data for the heterogeneous reactions (bl-b3) 
are derived from previous literature: combustion [37], steam 
gasification [38], carbon dioxide gasification [39]. 

Combustion of volatile products includes the reactions for 
primary tar and refractory tar (both modeled as hydrocarbon 
CH 1522 Oo.o 22 s [40]), with molecular weight equal to 94 [30,31], 
methane, carbon monoxide and hydrogen: 


(TO CH 1 . 522 O 0 . 022 s + 0.867O 2 ^ CO + 0.761H 2 O, 

CH 4 + 1.50 2 ^ CO + 2H 2 0 (cl -c2) 


2CO + 0 2 -32C0 2 , 2H 2 + 0 2 -52H 2 0 


(c3 —c4) 


(T 2 ) CHi. 522 O 0 . 0228 + 0.867O 2 ^CO + 0.761H 2 O (c5) 

The reaction rates are modeled as in Ref. [40], with the 
simplifications already introduced in Refs. [30,31] (first-order 
with respect to fuel (cl,c2,c5) and oxygen (c3)): 

r j = Aj exp ^ TCiCo 2 , j = cl, c2, c5; i = T 1 ,CH 4 ,T 2 (8a-8c) 

r c3 = A c3 exp f CcoCojC^o, r c4 = A c4 exp ^ c H2 c 02 

(9 -10) 

The corresponding kinetic constants are derived from 
Ref. [40] (cl,c2,c5), [41] (c3) and [30,31] (c4). The water gas shift 
reaction is also modeled: 

CO + H 2 0<—-C0 2 + H 2 , r wg = fe wg ^c co Ch 2 o-^—- j (wg, 11) 

with finite rate kinetics 

fe W g = A wg exp ^ ■ Kr = A e exp (42 —13) 

Different kinetic constants are used for the spatial zones of 
the reactor to modulate the catalytic action exerted by ashes, 
depending on their amounts [42]. Thus the kinetics reported in 
Ref. [43] is used for the emulsion phase and the splash zone, in 
the latter case corrected by a factor of 0.25 as already assumed 
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in Ref. [5]. The kinetics reported in Ref. [44], valid in the 
absence of significant catalysis, is used for the bubble phase 
and the freeboard (the equilibrium constant in all cases is 
derived from Ref. [45]). 


2.2. Hydrodynamics 


The hydrodynamic model is made by semi-empirical corre¬ 
lations for the variables of interest, such as the minimum 
fluidization velocity, u mf , [46], the bubble diameter, d b , [47] 
which also requires the bubble diameter at the top, d bm , [42], 
and bottom, d b0 , [46], of the bed, the single isolated bubble 
velocity, u br , [48], the volumetric fraction of the bubble phase, 
£ b , with the expansion factor B [42] and the bed height, H B , 
obtained from a mass balance on the solid phase: 


1.75 / diu mf p g \ 2 150(1 - r m f) ( d,u mfPg \ _ d ?Pg(ft ~ Pg)s 
4f«>i V J + 4f ?>i V M g y — 1 " 


^t>m db 
dbm — dbo 

d bm — min jo.652 





Iz=H b 



(15 -16) 


dbo — 2.78 


Umt ' 1 , u br = 0.711 y/gd^, 


£b = 1 


(17 -19) 


0738 „1.006 


B = 1 + 10.978 


(u 0 - u mf f ft 1 

<f 37 P ° g 126 

He 


W, = Sp,(l - E mf ) J (l-Eb)dz 


(20 - 21 ) 


The mass transfer between the emulsion and the bubble 
phase, driven by concentration gradients, is described ac¬ 
cording to Ref. [24] with the introduction of a global coeffi¬ 
cient, K be , consisting of the contributions of two series 
processes related to the exchange between bubble and cloud 
phase, K bc , and between cloud and emulsion phase, K ce : 


k qo. D° 5 g ° 2 

Kbe K bc + K ce ’ bC d b + ' d 225 


K,. = 6.77 


Dfmf^br i 


d 3 b 




(22 -24) 


The total mass flow rate is determined by the inlet gas 
stream, p g Su 0 , and the vapor/gas stream generated from wood 
drying and decomposition, F 0 (l - Y c ): 


PgSuo — p^Suo + F 0 (l — Y c ) (25) 

where Y c is the char mass fraction on a moist wood basis. 
Finally, the property equations are derived from Ref. [49] for 
the gas viscosity and from Ref. [50] for the species diffusivities. 

2.3. Particle size distribution 

Comprehensive models of single particle pyrolysis are avail¬ 
able including the effects of anisotropy, shrinkage and shape 
[51-55]. However, all the models of particle population bal¬ 
ances are based on one-dimensional spherical particles 
(shape effects can be taken into account by means of an 


equivalent diameter and a sphericity factor [34]). Moreover 
char particles preserve the same shape and size of the corre¬ 
sponding wood particles with a proportional decrease in the 
intrinsic density. These assumptions are also made in this 
study. 

The biomass particles present a fractional size distribution, 
p5, that should be assigned. More precisely, the lognormal 
function is selected [19]: 


Po = 


1 

d\/27rln a 


exp 


l/ln(d/ M )\ 2 ' 
2\ In a ) 


(26a) 


where fi is the mean and a is the standard deviation. The 
function is re-normalized over the range of particle sizes, 
df— d max , that is, the actual range of the particles fed, so that it 
shows a probability of 1: 


Po = - (26b) 

I Pod(d) 


Based on previous findings [17] and following the sche- 
matization used in other studies [18,20,21], it can be assumed 
that two contiguous size classes exist corresponding to par¬ 
ticles indicated as mothers and fines. A mother particle can be 
fragmented and can be consumed by chemical reactions and 
abrasion until reaching the upper size limit of the fines, d f , and 
becomes a fine itself. Thus the two particle size classes 
correspond to one range from 0 to d f for the fines and from d f to 
the maximum size, d max , of the mother class. The two classes 
present a continuous distribution of sizes within the respec¬ 
tive range. 

Once the char particles are formed, they undergo primary 
fragmentation that modifies the fractional particle size dis¬ 
tribution from po to Pq, although the shape remains 
unchanged: 

p'(d) = K f 1/3 p 0 (Kj /3 d) (27) 

where K f is a constant taken equal to 2 as previously done for 
the fuels of interest [20]. 

Both classes of char particles experience the heteroge¬ 
neous reactions of combustion and gasification where the 
density of the particle is constant but the size is shrinking 
with a rate R' [m/s]: 

R = +(+1 + 1+2 + + 3 ) (28a) 

The abrasion rate, R", of mother particles is described as in 
Ref. [56]: 

R" = y(ho-u mf ) (28b) 

where the parameter K a is equal to 3.5 x 10~ 7 , which is the 
range (2—7) x 10~ 7 reported for coal chars [18,20], leading to 
the formation of fines in accordance with an assigned frac¬ 
tional particle size distribution, p 0F , [17]: 


Pof — 


60/d 
df \d f 


Particles can be elutriated with a rate 


(29) 
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Rf^Ke,^, i = F, M (30) 

where P; is the particle size distribution (kg/m) of fines (i = F) or 
mother (i = M) particles and K el the elutriation constant. K el 
[57] essentially depends on the gas velocity, u 0 , and the ter¬ 
minal velocity, u t : 

K el = 0.011 ft /l-’^) (31) 


gasification (Eq. (28a)). The initial conditions again foresee a 
null distribution while the boundary condition requires the 
continuity between the two particle classes: 

P F (0, d) = 0. P F (t, d f ) = P M (t, d f ) (40 -41) 


2.4. Mass conservation equations for the emulsion 
phase 


where the terminal velocity u t is modeled as in Ref. [58]: 


u t = u t 


d = d 


1 -V3 


Mg I Pc Pg ) 9 


u t = 


18 0.591 


Pc Pg ) 9 


1/3 


(32 -34) 


The model equations for the perfectly mixed emulsion phase 
consist of the total mass conservation (42) with the reaction 
term (43) and the definition of the gas volume (44), the con¬ 
servation of single gaseous species (45a—h), the state equation 
(46), the definition of the gas density (47), and the mass of char 
particles for the mother and fine classes (48a-b) (the tem¬ 
perature is referred to the bed value): 


The mathematical model for the particle population con¬ 
sists of two unsteady one-dimensional (particle diameter) 
partial differential equations that are hyperbolic in character. 
The balance for the class of the mother particles is expressed 
in terms of the particle size distribution, P M , as: 


d Hb 

Ve dfP® = SUmf ( PAir “ p s) + zL / K ^,i f b(c[ - c^Mjdz 
1 0 

+ F 0 (l-Yc) + ]TX k -AF* 

k=l 


(42) 


® p _ p v n' P eI ® pp | foci 

c-P m - to icPo - k m - gd KF M H -d~ 

where the term on the left hand side is the accumulation. The 
contributions on the right hand side represent, in the order, 
the feeding rate, the elutriation rate, the convective term ac¬ 
counting for the displacement from larger to smaller sizes 
(owing to chemical reactions and abrasion) and the cumula¬ 
tive term for the chemical reaction and the abrasion rates. 
More precisely, R is the result of the contribution produced by 
chemical reactions (Eq. (28a)) and abrasion (Eq. (28b)): 

R = R' + R" (36) 

Initially the variable P M is assumed to be null over the 
entire range of the mother size class and the boundary con¬ 
dition states that there are no particles at the maximum size, 

dmax- 


d m ax df 

Rea = j r{j ^M c d(d) + I rj| ^M c d(d) (43) 

d f C 0 


Hb 

V e = Sf mf /(l-r b )dz 
0 


(44) 


d g FoYjo / o e\ ^et,k \r e 

Ve^rrC; — ^ + Su m f (q - Cj) + y '’ kr k 

1 k=l c k=4 


e dt 


H b 

+ J K be ,i£b(cf-c^dz-AFi, 

o 

i = 0 2 , CO, C0 2 , H 2 , H 2 0, CH 4 , T f , T 2 

(45 a —45h) 


Pm( 0, d) = 0, P M (t, d max ) = 0 (37 -38) 

The second partial differential equation is written for the 
class of fine particles in terms of the particle size distribution, 
P F , as: 

dmax 

f) r f) QD'P 

Ft Pr = Por J PuKd(d)-R?- M R'P f+ ^ (39) 

df 

The term on the left hand side is again the accumulation. 
At the right hand side, the various contributions, in the order, 
represent the production term resulting from the abrasion of 
the mother particles, the elutriation rate, the convective term 
taking into account the displacement from larger to smaller 
sizes, and the reaction rate. The production, due to mother 
particle abrasion, is expressed as the product of the total mass 
generated from abrasion and the resulting particle size dis¬ 
tribution which is assigned according to Eq. (29). The elutria¬ 
tion rate is again expressed by Eq. (30) while the reaction 
rate takes into account the contributions of combustion and 







(46 -47) 


(48a -48b) 


In addition to the mass exchange between emulsion and 
bubble phase, originated from the differences in species con¬ 
centrations and already taken into account through the global 
coefficient K be , the net flow, AF, [42] is also modeled. This is 
due to variations in the gaseous mass of the emulsion phase to 
guarantee a constant velocity, equal to the minimum fluid¬ 
ization velocity, following the generation of volatile species 
from wood drying and decomposition and is expressed as: 


AF: 



AF* > 0 


(49) 


The set of Eqs. (42,45a-45h) requires assigned initial 
conditions. 
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2.5. Mass conservation equations for the bubble phase 


The model equations for the one-dimensional, plug-flow 
bubble phase consist of the total mass conservation (50), the 
conservation of single gaseous species (51a-51h), the state Eq. 
(52) and the definition of the gas density (53) (the temperature 
is referred to the bed value): 


9 

at 



0Z ( £ bUbPg) + E Kbe.ifb (Cf 


c?)m,+ 


AF* 

SHb 


(50) 


f) f) 10 

- (EbC?) = (fbUbC?) + J2 "U^b + K be ,i£ b (cf 

k=4 

i = 0 2 ,C0,C0 2 ,H 2 ,H 2 0,CH4,Ti,T2 


«f) + 


AF ; 

SHb’ 


(51a —51h) 


c n 2 =S-E c ?’ Pg = E c ? M ' ( 52 ~ 53 ) 

i i 

with assigned initial conditions. The boundary conditions 
correspond to known values of the inlet mass fluxes (at z = 0). 


2.6. Mass conservation equations for the splash and 
freeboard zones 

The model equations for the one-dimensional, plug-flow 
splash and freeboard zones, with variable cross section S S f(z), 
consist of the total mass conservation (54), the conservation of 
single gaseous species (55a-55h), the state equation (56), the 
definition of the gas density (57) and the mass conservation 
equation for char (58) with average size of the char particles 
(59) (the temperature is referred to the bed or the freeboard 
value): 


3. Analytical solutions of the particle 
population balances 

To develop an analytical solution of the model equations for 
particle size distribution, in the first approximation, it is 
assumed that the conditions of the surrounding gas (velocity, 
temperature and species concentration) are assigned at con¬ 
stant values, so that reaction, abrasion and elutriation rates 
are also constant. The analytical solution is developed with 
the same approach detailed in Ref. 19], although significant 
differences exist between the two mathematical models. 
Therefore modifications are introduced in the solution pro¬ 
cedure in order to take correctly into account the new fea¬ 
tures. In this study two size classes of particles are considered 
versus the single class used in Ref. [19], Hence the solutions 
should be determined for two partial differential equations 
(Eqs. (35) and (39)), which are coupled through the time- 
dependent boundary condition (Eq. (41)), instead of the 
simpler constant value at the boundary imposed in the pre¬ 
vious work. Moreover, the present treatment introduces some 
additional features: a) a more realistic description of the het¬ 
erogeneous reaction rates is used (Eqs. (2) and (3)) instead of a 
the simple expression A'/(d + B'), where A' and B' are assigned 
constants, b) the effects of attrition are taken into account and 
c) a more general treatment of the elutriation process is pro¬ 
posed instead of a simple on/off process where all fuel parti¬ 
cles smaller than a critical size are elutriated whereas the 
larger particles remain in the system. 

The following expressions are obtained for P M and P F : 

P M (d,t) = -F 0 Y c |- j ^-u^t + yid(d)jexp^^||d(d)jd(d') 

(60) 





( B4 ) P F (d, t) = P F1 (d, t) + P F2 (d, t) 


(61) 


s(<f) = |(Wf) + f f + f>rf, 

-bSF OZ *5 d S F Pc j—J 

i = 0 2 ,C0,C0 2 ,H 2 ,H 2 0,CH 4 ,Ti,T 2 


€=g-E c ? F > />f = E c ? FM i 


(55a —55h) 
(56 -57) 


P F1 (d,t) = -^J P fu(t+ J ^ d (d)]rnJt + j ^d(d) 


xexp 


7 d(d) j d(d') 

V* ■ ■ J 


K d S 

WjR' 


(62a) 


ofi ~ “s^az (SsfUsfP * F) - E 


k=l U SF 


dsF Pc 


(B8) PF2 (d,t) = -(|) ^P M (d f ,t)exp 


( d f \ 

v° / 


(62b) 



(59) 


Boundary conditions, at the bottom of each zone, require 
the continuity of the mass fluxes taking into account that char 
particle transport occurs via elutriation. 


where ifi F (t) represents the mass flux of fines produced by 
abrasion of mother particles according to the following 
equation (Eq. (63)): 

dmax 

m F (t)= j P M RSd(d) (63) 

i, 

Eqs. (60)-(62) provide the following steady-state distribu¬ 
tions for t —» oo: 
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dmax / H \ 

,(d) = -F 0 Y c | j JexpM^|d(d)jd(d') 


P F(d) = Pn(d ) + PF2(d) 


Pn(d) = -R|) / /f eXP (/If d(d) ) ‘ 

d, d \ d / 


(64) 

(65) 
| d(d') (66a) 


Pf2 (d) 



3 R'(d f ) 

R'(d) 


P M (d f )exp 



(66b) 


In the problem schematization it is assumed that the 
mother and fine classes of particles are contiguous. However, 
it has been demonstrated [17,18,20,21,23] that in fluidized bed 
reactors attrition essentially occurs by means of an abrasive- 
type mechanism and thus, between the son particles (a sub¬ 
class of the mother particles, produced from the fragmenta¬ 
tion and abrasion processes) and the fines produced, the range 
of intermediate-sized particles is almost empty in terms of 
mass. This circumstance provides the basis for an analysis of 
the characteristic evolution times for the two size classes 
introduced which, given the large size differences, are ex¬ 
pected to show highly different values. Indeed, as discussed in 
the following, it is found that, for the conditions typical of 
fluidized-bed gasification, the characteristic times for the 
fines are several order of magnitude shorter than those of the 
mothers. Consequently it is reasonable to assume that the 
class of fine particles exhibits a quasi-steady behavior with 
respect to the class of mother particles. In this way, the un¬ 
steady equation for the population balance for the mother 
class (Eq. (35)) can be coupled, through the boundary condi¬ 
tions (Eq. (41)), to the steady-state version of the population 
balance equations of the fines (Eq. (39)) and the solution given 
by Eq. (60) with that given by Eq. (65). The validity of the quasi¬ 
steady approximation will be discussed in relation to the 
range of typical operating conditions of fluidized-bed gasifers. 


4. Numerical solution 

The assumption that the environment external to the particles 
is at uniform and constant conditions in terms of temperature, 
species concentrations and velocity is no longer valid when 
the particle population equations are coupled with chemical 
reactions and the heat and mass transfer processes for the 
fluidized-bed gasifier. On the other hand, given the complexity 
of the conservation equations for the latter, a numerical so¬ 
lution should be anyway implemented. A semi-implicit finite- 
difference procedure is applied here. For each time step, the 
numerical solution for the population balance equations is 
carried out first, using the last available values for the variables 
obtained from the equations for the different zones of the 
fluidized bed. Then these are solved using the solution already 
determined for the particle size distribution. The two stages of 
the solution procedure are repeated for successive time in¬ 
tervals until steady-state conditions are attained. 


The numerical solution of the population balance equa¬ 
tions is based on the assumption of quasi-steady behavior, 
that is, the analytical solution (65) of the steady formulation of 
the equation for the fine particle class is combined with the 
numerical solution of the unsteady equation (35) for the 
mother particles. This is based on the operator splitting pro¬ 
cedure [59] and finite differences approximations of the con¬ 
vection terms using the Flux-Corrected-Transport (FCT) 
scheme with reference to the explicit SHASTA phoenical 
version [60], The procedure foresees two stages corresponding 
to convective transport and the remaining contributions at 
the right hand side of the equation. Thus, for each time step, 
after convective transport, the odel5s integrator of the Matlab 
platform is applied. Then, using the last available values for 
the abrasion rate and the boundary condition (41), the 
analytical solution (65) is applied for the class of fine particles. 
In this stage of the computational procedure, the conditions 
external to the particles are given by the last available values 
of the solution of the equations for the fluidized-bed reactor. 

As for the different spatial zones of the fluidized-bed 
reactor, the convective terms of the equations for the bub¬ 
ble, the splash and the freeboard zone are discretized by 
means of the finite difference backward scheme. In this way, 
each partial differential equation is approximated by a system 
of N ordinary differential equations (N is the number of spatial 
cells). A staggered grid is employed where the process vari¬ 
ables are positioned at the center of the cell and the velocity at 
the cell boundaries. The resulting system of ordinary differ¬ 
ential equations together with the ordinary differential 
equations for the emulsion phase and all the algebraic equa¬ 
tions are solved using again the odel5 integrator of the plat¬ 
form Matlab. In this stage of the computational procedure, the 
particle size distribution coincides with the last available. 


5. Results 

Wood gasification with air in a fluidized bed is studied. The 
results consist of two main parts. In the first one, the dy¬ 
namics of particle size distribution are discussed. In the first 
place, the analytical solutions of the particle population bal¬ 
ances are applied to simulate the effects of the composition of 
the surrounding gas and the feed size. In this way, the con¬ 
ditions which permit the application of the quasi-steady 
approximation of the fine with respect to the mother class 
are also determined. In the second part, results are discussed 
of the simulations of the coupled numerical solution of the 
particle population balance model with the transport equa¬ 
tions for a bubbling fluidized-bed gasifier, providing a com¬ 
parison between predictions and measurements of the 
composition of the gas produced at steady-state conditions. 

5.1. Dynamics of particle population 

In this section results are presented about the dynamics of 
particle size distribution, assuming that the environment 
external to the particle is at assigned conditions. Typical 
operating conditions correspond to a fuel feed rate of 
0.28 kg s _1 with an air (kg) to fuel (kg) ratio of 1.8, a bed tem¬ 
perature of 1073 K and a fluidizing gas velocity of 1 m s \ The 
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bed is assumed to consist of silica sand particles (size 500 (im) 
with a sphericity factor of 0.86. For the reference case, the 
particle size distribution of the feed is assumed to be 
comprised (99.5%wt) in the range d min -d max with values of 10 
and 30 mm, respectively, whereas the limit diameter, d f , 
positioned at the boundary between the classes of mother and 
fine particles, corresponds to 2.4 x 1CT 3 mm [20]. The density 
of the char particles is 125 kg/m 3 . Two cases are considered for 
the chemical composition of the gaseous environment sur¬ 
rounding the particles: A) it is a mixture obtained from the 
pyrolysis gas and the inlet air (molar composition: X 0 2 = 0.141, 
X N2 = 0.5295, X co = 0.091, X C02 = 0.023, X H2 = 0, X H20 = 0.154, 
Xch 4 = 0.052, X T1 = 0, X T2 = 0.0092) or B) it is the typical gas 
composition at the exit of a bubbling fluidized-bed gasifier 
(molar composition: X 02 = 0, X N2 = 0.36, X C o = 0.14, 
X C02 — 0.12, X H2 = 0.10, Xpj2o 0.28, Xc^4 — 0.0, X FF = 0, X F2 — 0). 
These two atmospheres are indicated in the following as A 
and B. 

Three cases of mass transfer between gas and particles 
have been identified for fluidized-bed reactors depending on 
the relative size of the active particles and the bed particles 
[61]. For particles much smaller than the bed particles it is 
appropriate to use the Frossling correlation [62]. For particle 
sizes comparable to the bed particles, the correlation pro¬ 
posed by La Nauze and Young [63] is widely used which, for 
larger particles, is substituted by the correlation proposed by 
Prins et al. [64], Given the range of particle sizes typical of 
fluidized bed reactors (2-5 mm [65]), the chemical times 
evaluated with the above correlations for the problem under 
study show values coincident or only slightly shorter than 
those obtained by means of the Ranz-Marshall correlation. As 
already observed, in practical cases, the lower-size part of the 
mother class is almost empty, which also enhances the actual 
differences between the characteristic times for the two 
classes. Based on these considerations, it is important to 
establish the validity limits of the assumption of a quasi¬ 
steady behavior of the fine versus the mother class. 

To analyze the dynamics of the particle size distribution, 
the unsteady analytical solutions (60,61) are applied for the 
reference case. A grid of 1000 cells, for each of the two particle 
classes, is used for the numerical evaluation of the integrals (a 
grid of 100 cells for the fine class introduces deviations on 
output variables below 3%). The variables P M and P F , as 
computed for the atmosphere A, are shown for several times 
in Fig. 1A and B. As expected, the examination of the various 
global rates (the contributions in the conservation equations 
integrated over the particle size) confirms that for the mother 
class the dominant process is the chemical reaction rate that, 
at steady conditions, practically balances the inlet rate of the 
fragmented char particles. The contributions due to abrasion, 
convective outflow and elutriation rates are smaller by 
approximately two, eight and six orders of magnitude, 
respectively. The main inlet contribution for the fine class is 
due to the abrasion rate of mother particles which, at steady 
conditions, is again balanced by the chemical reaction rate 
(the convective inlet rate and the elutriation rate are smaller 
by about six orders of magnitude). Both the mass, for a given 
particle size, and the different global rates continuously in¬ 
crease in time until steady conditions are achieved for times 
around 260s (corresponding to the attainment of 95% of the 



d x 10 3 [mm] 

Fig. 1 — Particle size distributions for the mother class, P M , 
(A) and the fine class, P F , (B) for several times and 
atmosphere A (reference case). 


total final char bed loadings for the two particle classes). Given 
the coupling resulting from the convective outflow and the 
abrasion rate from the mother class to the fine class, the times 
needed to achieve steady values of the related char bed 
loading are the same. 

The P M and P F curves provide information about the mass 
of char particles in the bed of a given size. Thus, in general, 
both the mass and the size range may change in time. Indeed, 
Fig. 1A shows that, while the mass increases with time, the 
particle size range for the mother class tends to be displaced 
towards successively smaller values. This can be better seen 
from Fig. 2, which shows the corresponding fractional particle 
size distribution, p M (that is, P M referred to the total mass), at 
steady conditions and the corresponding values for the char 
feed, p 0 , and the fragmented char particles, p 0 '. For the mother 
class the times needed to achieve the steady conditions, in 
terms of mass and particle size range, are the same and 
approximately correspond to the characteristic chemical 
times that, for the sizes of interest, are of the order of hun¬ 
dreds of seconds. Instead, the P F curves reported in Fig. IB 
show that, while the mass continuously increase in time, the 
size range for the fine class is invariant already at short times. 
Early process dynamics show that a steady range of particle 
size for the fine class is attained within times of the order of 
milliseconds. Again, these times are essentially those of the 
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Fig. 2 - Steady-state fractional particle size distributions 
for the mother class, p M , and the fine class, p F , for the 
atmosphere A (solid line) and atmosphere B (dashed line) 
and the reference case. For comparison purposes, the 
fractional particle size distribution of the feed, p 0 , and the 
fragmented particles, p 0 ', and the abraded particles, p 0F , are 
also included. 


chemical reactions. Also, in this case there is a displacement 
of the particle size range towards lower values as shown by a 
comparison between the fractional mass size distribution at 
steady conditions, p F , and the corresponding value for the fine 
distribution from abrasion, p 0F in Fig. 2. In other words, 
although inlet and outlet global mass feed rates change over 
long times, as determined by the characteristic times of the 
mother particles, given that the size range of inlet particle is 
constant in time, a steady distribution is achieved within the 
times required by the chemical reactions to occur. These, 
being the particle sizes very small, are much shorter than 
those typically associated with the size evolution of the 
mother class. 

The process dynamics remain qualitatively the same for 
the atmosphere B and the relative importance of the various 
contributions in the conservation equation is also left unal¬ 
tered. The steady fractional mass distributions reported in 
Fig. 2 indicate that the particle sizes of the mother class are 
barely displaced towards smaller values whereas for the fines 
there are no appreciable changes. The curves of the char bed 



Fig. 3 — Char bed loading versus time for the mother class, 
W M , and the fine class, W F , for the atmosphere A (solid line) 
and atmosphere B (dashed line) and the reference case. 


loading versus time (Fig. 3) show that higher values are 
reached for the atmosphere B (factors of about 1.5 and 8.5 for 
the mother and fine classes) and that the times needed to 
achieve steady conditions are also longer (about 360s versus 
260s). These features can be explained by observing that the 
fastest chemical reaction rate for the atmosphere B (steam 
gasification) presents characteristic times that, for a kinetic or 
a mixed kinetic-diffusive control, are significantly longer than 
those of the oxidation rates of atmosphere A. The qualitative 
features of the particle size distribution remain the same as 
the gas velocity and the temperature of the surrounding 
environment are varied (the abrasion and elutriation rate in¬ 
crease as the velocity is increased or the temperature is 
decreased). 

To further investigate the dynamics of particle size distri¬ 
bution, feed sizes are varied in the ranges d min —d max (99.5%wt) 
of 25—50, 20—40, 5—20,1-5 and 0.1-1 mm. A grid consisting of 
1000 cells, for each of the particle class, gives accurate 
computation of the integrals. The qualitative features of the 
process again remain the same as discussed for the reference 
case (d min —d max = 10-30 mm). As expected, the time needed to 
attain steady conditions and the steady char bed loadings 
show successively lower values as the feed sizes are progres¬ 
sively reduced. The increase in the rate of the heterogeneous 
reactions of char conversion, the dominant contribution in the 
particle size population balance, is responsible for such a 
behavior. Also, the qualitative effects of the reaction envi¬ 
ronment remain the same as already reported for the refer¬ 
ence case. For instance, the time needed to attain steady 
conditions varies from about 720 to 1 s (atmosphere A) or from 
about 880 to 6 s (atmosphere B). The steady char bed loading 
for mother particles is also different with values in the range 
24-0.04 kg (atmosphere A) or 29-0.18 kg (B). The times needed 
to achieve steady conditions for the class of fine particles are 
the same as for the mother class (the steady char bed loading 
values also decrease although values are always very small). 

Examples of the steady profiles of the fractional particle 
size distribution for the mother class are plotted in Fig. 4 



Fig. 4 - Steady-state fractional particle size distributions 
for the mother class, p M , for the atmosphere A and various 
ranges d m i n -d max of the feed particles (solid line). For 
comparison purposes, the fractional particle size 
distributions of the fragmented particles, p 0 ', are also 
included (dashed line). 
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where, for comparison purposes, the initial distribution for 
the fragmented particles is also shown. The two profiles are 
similar although, due to chemical reaction and abrasion ac¬ 
tivity, the percentage of smaller particles for the former in¬ 
creases significantly especially when larger feed particles are 
considered. Indeed, the maximum p M is always located in 
correspondence of the descending zone of the left-hand side 
of the p 0 ' distribution and attain a lower value. Obviously, the 
dominant term, at steady conditions practically balancing the 
feed rate, is always due to the chemical reaction rate. How¬ 
ever, by decreasing the feed size range, at steady conditions, 
the global elutriation rate and the convection rate increase by 
about four and two orders of magnitude (atmosphere A) or five 
and two orders of magnitude (atmosphere B) (information on 
the steady profile of the elutriation rate distribution can be 
obtained from Fig. 5, where maximum values are increased up 
to factors of about 7 (atmosphere A) or 28 (atmosphere B)).The 
steady abrasion rate decreases by approximately one order of 
magnitude. It is important observing that, given the constant 
particle size distribution from the abrasion rate p 0F , the steady 
fractional particle size distribution for the fine particles re¬ 
mains unvaried with respect to the profile reported in Fig. 2. 
Indeed the chemical reaction rate again practically balances 
the abrasion rate for all the feed size ranges investigated. 
However, as the feed particles are displaced towards smaller 
values, given the constant distribution in the size of fine 
particles and the progressive diminution in the char bed 
loading, the elutriation rate decreases, remaining from about 
six to five orders of magnitude slower (see Fig. 5 for the steady 
profile of the elutriation rate distribution which shows re¬ 
ductions up to factors of about 9 (atmosphere A) or 3 (atmo¬ 
sphere B)). The steady convective transport rate are also still 
slower by about six or three orders of magnitude. 

Computations of the dynamics of particle size distribution 
for all the cases examined, which approximately span the 
typical conditions established in fluidized-bed gasifiers, pro¬ 
vide the same results by using either the fully unsteady solu¬ 
tions or the quasi-steady-solutions, indicating that the latter 
can be safely applied. To further investigate the applicability of 



Fig. 5 - Steady distribution of the elutriation rate for the 
atmosphere A (solid line) and atmosphere B (dashed line) 
and various ranges d min -d max of the feed particles: the 
elutriation rate increases and decreases, according to the 
direction of the main arrows, for the fine and mother 
classes, respectively. 


this treatment the case is considered of a very small feed size 
range corresponding to d min —d max = 0.03—0.1 mm, that is, for 
conditions of a very close coupling between the two size 
classes (grids with 1000 cells again provide an accurate nu¬ 
merical evaluation of the integrals). Fig. 6A and B summarizes 
the results in terms of particle size distribution (at several 
times) and char bed loading profiles for the two classes of 
particles and the atmosphere A. The steady distribution of the 
particle size for the mother class is always dictated by a bal¬ 
ance between the feed rate and the chemical reaction rate with 
the elutriation, abrasion and convection rate smaller by about 
four, three and four orders of magnitude respectively. Given 
the small particle sizes and the fast rate of heterogeneous 
conversion, steady conditions are attained within very short 
times (about 0.1 s). Times needed to achieve steady conditions 
are longer (by a factor of approximately 2) for the class of fine 
particles. In this case, the reaction rate is still the most 
important outlet contribution (the steady elutriation rate is 
smaller by about six orders of magnitude). However, the inlet 
rate is represented not only by the abrasion rate of the mother 
particles but also by the convection term (boundary condition) 
which, at steady conditions, is smaller only by a factor of about 
2.7. In other words, the convection term is negligible for the 
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Fig. 6 - Fractional particle size distribution for the mother 
class, p M , (A) and fine class, p F , (B) and corresponding char 
bed loadings, W M and W F , versus time for several times 
and atmosphere A (feed range d min -d max equal to 
0.03 mm-0.1 mm). 
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Table 1 - Predicted (S) and measured (E) gas composition, 
six tests reported in Ref. [66]. 

as percent on a molar basis, and LHV and tar content of the gas for 

Test number 


1 


2 

3 



4 

5 



6 

Folgs- 1 ] 


0.16 

0.153 

0.108 


0.19 

0.182 

0.188 

Moisture content [%wt] 

23.5 


21.0 


23.0 


22.0 


25.0 


19.0 


ER 

0.32 


0.37 


0.47 


0.26 


0.36 


0.32 


Tb[K] 

1073 


1073 


1083 


1073 


1063 


1073 


T f [K] 

813 


823 


773 


873 


833 


803 


Carbon conversion [%wt] 

84.4 


91.2 


71.5 


76.4 


87.3 


92.3 



E 

S 

E 

S 

E 

S 

E 

S 

E 

S 

E 

S 

CO 

14.0 

14.4 

13.0 

15.2 

10.0 

10.2 

13.0 

15.8 

13.0 

15.4 

18.0 

17.8 

co 2 

13.5 

14.1 

15.0 

13.6 

12.0 

14.3 

15.0 

13.9 

15.0 

13.2 

13.5 

12.3 

ch 4 + c 2 h 4 

4.2 

4.7 

4.3 

4.6 

3.5 

2.1 

4.3 

5.6 

4.3 

4.2 

6.8 

4.7 

h 2 

7.0 

10.0 

9.5 

10.4 

8.0 

4.9 

9.5 

12.0 

9.5 

10.0 

9.5 

11.8 

n 2 

61.3 

69.0 

58.3 

56.1 

66.5 

68.4 

58.3 

52.7 

58.3 

57.2 

52.2 

53.3 

h 2 o 

8.0 

26.1 

7.7 

24.1 

27.0 

26.4 

21.6 

26.5 

24.0 

26.6 

0.24 

23.1 

TAR [g m~ 3 ] 

3.73 

6.99 

7.16 

7.09 

2.99 

3.17 

2.01 

8.32 

2.01 

6.34 

9.98 

7.36 

LHV [MJ m~ 3 ] 

4.3 

4.6 

4.6 

4.7 

3.7 

2.6 

4.6 

5.3 

4.6 

4.5 

6.3 

5.2 


population balance of the mother particles but it is important 
for the fine particles. It increases in time and causes a pro¬ 
gressive variation in the shape of the size distribution of fine 
particles. At very short times this reproduces the same quali¬ 
tative trends shown in Fig. IB. However, as the convection rate 
increases, the mass of particles with sizes displaced towards 
larger values becomes progressively higher. The time profile of 
the char bed loading clearly shows that there is a first zone, 
approximately corresponding to times needed by mother 
particles to attain steady conditions, essentially determined 
by the abrasion rate of mother particles and a second period 
determined by both abrasion and convection rates of mother 
particles. Similar results are also obtained for the atmosphere 
B although the time needed to achieve a steady state is longer. 
For this range of small-sized feed particles the quasi-steady 
solution over-predicts the particle size distribution of fine 
class with deviations that are time and size dependent. 
Maximum deviations, corresponding to about 40%, are 
observed around the boundary between the fine and the 
mother class for a time approximately coincident with that of 
the steady solution of the mother class. Deviations for the 
smaller sizes are approximately around 10%. In all cases, they 
tend to decrease as steady conditions for the class of fine 
particles are approached. Therefore, in the case of a close 
coupling, the assumption of a quasi-steady behavior of the fine 
with respect to the mother class cannot be applied. 

The analytical solution of the particle population balance 
has also been applied to evaluate the accuracy of the nu¬ 
merical solution using the reference case. Simulations have 
been made using 100, 200, 300, 400 and 500 cells with time 
steps of 0.01 s for Eq. (35). Deviations between the numerical 
and the analytical solutions, in terms of particle size distri¬ 
bution and char bed loading, at very short times (below 10 s) 
reach values up to 20-30% that appear to be insensitive to the 
number of cells. Then deviations on the particle size distri¬ 
bution decrease to values below 1% already for a grid of 300 
cells (about 10% for a 100 cell grid). A similar trend is observed 
for the deviations in the char bed loading which also shows 
lower values. Thus, apart from the very initial transients, the 
numerical solution gives accurate predictions already with a 
relatively coarse grid. 


5.2. Bubbling fluidized bed simulations 

The coupled model for the particle size distribution and the 
bubbling fluidized bed is applied to simulate the laboratory- 
scale system developed in Ref. [66]. It is made by a bottom 
vessel with a height of 0.6 m and diameter of 0.06 m, where the 
expanded bed and the splash zone are located, and un upper 
vessel with a height of 0.2 m which also present an increase in 
the diameter (over a length of 0.05 m) from 0.06 m to 0.12 m, 
which represents the freeboard zone. The bed material is silica 
sand (0.5 kg) with sizes between 320 pm and 500 pm. The wood 
feed rates are approximately in the range 0.1 g s~ 3 —0.2 g s~ a 
with equivalent ratios (ER) between 0.26 and 0.47 (corre¬ 
sponding air (kg) to fuel (kg) ratios of 1.4-2.4), the bed tem¬ 
peratures, T b , around 1063-1083 K and freeboard 
temperatures, T F , around 773-873 K. Pine wood particles sizes 
are between 0.8 mm and 4 mm and the initial density is 
450 kg m~ 3 which, in accordance with the stoichiometry of the 
pyrolysis reaction, gives rise to a char density of 96.75 kg mr 3 . 
The maximum diameter for the class of fine particles, d f , is 
again assumed equal to 2.4 x 10~ 3 mm [20]. The carbon 



Fig. 7 - Measurements (symbols [66]) and predictions 
(lines) of the dry gas composition, on a molar basis, versus 
the equivalent ratio, ER (wood feed rate 0.12 g s -1 , 
moisture content 21% wt, carbon conversion 84% wt, 

T b = 1073 K, T F = 823 K). 
















134 


BIOMASS AND BIOENERGY 62 (2OI4) I 2 3 —I 3 7 



Fig. 8 - Measurements (symbols [66]) and predictions 
(lines) of the LHV of the gas and the tar content of the gas 
versus the equivalent ratio, ER (the measured LHV values 
include those given in Ref. [66] (open symbols) and those 
evaluated from the gas composition (standard conditions) 
reported in the same study (full symbols)) (conditions as in 
Fig. 7). 


conversion, evaluated from the elemental composition of 
wood and the composition of the gas provided for six typical 
gasification experiments is between 70 and 90%. In the simu¬ 
lations it is assumed that the unconverted carbon belongs to 
the char product with a corresponding decrease in the feed 
rate. The initial conditions foresee that the gas phase consists 
of air at the assigned temperature for the bed and the free¬ 
board. Then wood is fed to the gasifier and the conversion 
process begins. Feeding takes place in the emulsion phase and 
the excess gas/tar produced, with respect to the conditions of 
minimum fluidization, is uniformly distributed across the 
bubble phase. 

Similar to previous models of fluidized-bed gasifiers [4-12], 
the experimental validation of the model is made considering 
the composition of the gas, as this is largely dependent on the 
chemical kinetics and transport phenomena [67], Indeed, for 
these gasifiers, the temperature is approximately uniform but 
rather low. Hence heterogeneous chemical equilibria will 


usually not be attained and only serve to set the trend of the 
variations in gas composition with the operating parameters. 
Homogeneous equilibrium is more likely to be established, as 
the water gas shift reaction is essentially in equilibrium above 
1100 K. 

An acceptable agreement between measured and pre¬ 
dicted properties of the gas for six typical experiments re¬ 
ported in Ref. [66] is obtained as shown by the data listed in 
Table 1. It is plausible that there is a misprint in the reported 
experimental values for the H 2 0 yields for the experiments 

N. l and N.2, as these values are lower by a factor of about 
three with respect to those of the other experiments and the 
predictions of the model. 

A further comparison between predictions and measure¬ 
ments can be made through Fig. 7 (gas composition, as a 
percent on a molar dry basis) and Fig. 8 (gas lower heating value 
(LHV) and tar content of the gas) for ER values in the range 

O. 25-0.45 (wood feed rate 0.12 g s _1 , moisture content 21%, 
carbon conversion 84%, T B = 1073 K, T F = 823 K). As ER is 
increased, the CO, C0 2 , CH 4 and H 2 concentrations decrease as 
a consequence of the successively higher amounts of nitrogen 
which prevails upon the effects deriving from an enhanced 
activity of the combustion reactions which are responsible for 
the reduction in the tar content of the gas. In general, the model 
provides good quantitative results although the predictions 
indicate a stronger dependence of the H 2 concentration on ER. 
The maximum steady elutriation rate distribution increases 
with ER with values in the range (1-2) x 1CT 5 kg nT 1 s _1 for the 
mother class and corresponding values for the fines lower by 5 
orders of magnitude. 

Acceptable agreement between predictions and measure¬ 
ments of the gas composition is also obtained for various bed 
temperatures as shown in Fig. 9 for the gas composition and 
Fig. 10 for the gas LHV and the tar concentration (wood feed 
rate 0.67 kg/h, moisture content 21%, carbon conversion 84%, 
ER = 0.3, T F = 923 K). As the bed temperature is increased, the 
steam gasification of char becomes more favored with respect 
to C0 2 gasification (higher activation energy) with an increase 
in the production of H 2 . Moreover, the forward path of the 



Fig. 9 - Measurements (symbols [66]) and predictions 
(lines) of the dry gas composition, on a molar basis, versus 
the bed temperature, T B (wood feed rate 0.19 g s' 1 , 
moisture content 21% wt, carbon conversion 84% wt, 

ER = 0.3, T f = 823 K). 
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Fig. 10 — Measurements (symbols [66]) and predictions 
(lines) of the gas LHV and tar content of the gas versus the 
bed temperature, T B (the measured LHV values include 
those given in Ref. [66] (open symbols) and those evaluated 
from the gas composition (standard conditions) reported in 
the same study (full symbols)) (conditions as in Fig. 9). 
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water gas shift reaction becomes less favored, with an in¬ 
crease in the production of CO (and a diminution in the pro¬ 
duction of C0 2 ). The CH 4 yields remain approximately 
constant. The model over-predicts the yields of H 2 . However, 
as observed in Ref. [66], the measured yield of this products 
vary over a rather wide range that, for the temperatures of 
interest here, is about 8-14% and is comparable with the 
values given by the model. Finally, although the predicted tar 
content of the gas is close to the measured values, the model 
does not predict the slight decrease as the temperature is 
increased. This can be attributed to steam reforming that is 
not included in the reaction network considered here. Again, 
the elutriation rate distribution decreases as T B is increased 
with maximum steady values diminished by one order of 
magnitude over the range of interest. 


6. Conclusions 

A comprehensive unsteady model for the distribution of par¬ 
ticle sizes in fluidized-bed reactors has been developed 
including fragmentation, abrasion, elutriation and the chem¬ 
ical reactions of wood gasification. The particle distribution is 
divided into two size ranges: the class of fine particles 
contiguous to the class of mother particles. Simulations, car¬ 
ried out by varying temperature, composition and velocity of 
the surrounding gas and feed size range over conditions that 
span those typically established in fluidized-bed gasification, 
have been applied to estimate the relative contribution in the 
balance equations, in particular the elutriation losses. These 
become successively more important as the gas velocity is 
increased, the temperature is decreased and/or the feed size is 
decreased. An important finding of the simulations is that, for 
the conditions of practical interest, a quasi-steady approxi¬ 
mation can be made for the behavior of fine particles with 
respect to the class of mother particles. 

A mixed numerical (mother class) and analytical (fine 
class) solution of the quasi-steady population balance equa¬ 
tions is incorporated in a model for a bubbling fluidized bed 
biomass gasifier schematized according to the two-phase 
theory of fluidization, where the emulsion phase is assumed 
to be perfectly mixed and a plug-flow approximation is used 
for the bubble phase and the splash and freeboard zones. The 
coupled particle population balance and the fluidized bed 
models are used to simulate gasification of wood with air at a 
laboratory scale. Simulations made for the typical conditions 
of the experiments and by varying the equivalent ratio and the 
bed temperature show good agreement with the measure¬ 
ments for the gas composition (and LHV) and the tar content 
of the gas. 

Although the model presented in this study represents a 
significant advancement with respect to the current state of 
the art, there are still numerous aspects that deserve further 
developments. These include the formulation of more accu¬ 
rate sub-models for the processes of particle size reduction 
(primary and secondary fragmentation, abrasion) and the 
experimental determination of the limit particle size between 
mother and fine classes. Extensive experimental investiga¬ 
tion, which can provide detailed information on the particle 
size distribution on dependence on the operating conditions, 


is also highly desirable not only for an improved under¬ 
standing of the process but also for the experimental valida¬ 
tion of the model proposed. The mathematical model should 
be improved with a more accurate description of the reaction 
rates with the assumption of regimes that should vary, for 
instance, on dependence of the combustion or the gasification 
process. Also, effects associated with non-isothermal 
behavior may be important and should be properly taken 
into account. 


Nomenclature 

A pre-exponential factor [s -1 ] 

c species concentration [kmol m -3 ] 

d particle diameter [m] 

d f boundary diameter between the class of fine and 

mother particles [m] 

d max maximum particle diameter [m] 

dmin minimum diameter for 99.5% of feed particles [m] 

D diffusivity [m 2 s -1 ] 

D bed diameter [m] 

E activation energy [J kmol -1 ] 

ER equivalent ratio 

F 0 wood feed rate [kg s -1 ] 

g acceleration of gravity [m s -2 ] 

K a abrasion constant 

K f fragmentation constant 

H b expanded bed height [m] 

LHV gas low heating value [J m -3 ] 

M molecular weight [kg kmol -1 ] 

P particle size distribution [kg itT 1 ] 

R g Universal gas constant [J kmol -1 K -1 ] 

S area of reactor cross section [m 2 ] 

T temperature [K] 

t time [s] 

u velocity [m s -1 ] 

z spatial coordinate [m] 

W] sand mass [kg] 

X molar fraction 

Y mass fraction 

Greek symbols 
e void fraction 

,u mean value of the input fractional size distribution 

[m] 

,u g viscosity [Pa s] 

v stoichiometric coefficients for biomass primary 

pyrolysis 

v* stoichiometric coefficients for thermal cracking of 

primary tar 

p density [kg m -3 ] 

a standard deviation of the input fractional size 

distribution 

ip sand sphericity 

n pressure [Pa] 

Superscripts 

e Emulsion phase 

0 Input conditions 
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b Bubble phase 

SF Splash and freeboard 

Subscripts 

bl char oxidation 

b2 char gasification (C0 2 ) 

b3 char gasification (H 2 0) 

cl primary tar combustion 

c2 CH 4 combustion 

c3 CO combustion 

c4 H 2 combustion 

c5 refractory tar combustion 

c char 

E equilibrium 

F fine particles 

g gas phase 

M mother particles 

mf minimum fluidization conditions 

S solid phase 

SF splash and freeboard zone 

wg water gas shift 
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